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A dense phase of GaAs wires forms in the early stages of 
strained growth on GaP,assembling from elongated Stranski- 
Krastanow islands. The electron diffraction during growth is 
consistent with long, faceted GaAs islands that are anisotrop- 
ically deformed without dislocations. The lateral wire pe- 
riod and long shapes are not predicted by published models, 
though we conclude that the island orientation is picked out 
by facet energy inequivalencies not present in the analogous 
system of Ge islands on Si. 

PACS numbers: 81.15.Hi, 68.35. Rh, 68.35. Bs, 61.14.Hg 



Crystal growth on a substrate with a smaller lat- 
tice spacing can initially give small regular islands with- 
out crystal dislocations. Two examples are the tens-of- 
monolayers-high Ge islands that form on Si |] and InAs 
islands that form on GaAs pH. The layer-substrate 
system lowers its energy when the islands relax towards 
their larger lattice constant at the expense of both locally 
straining the substrate and increasing the surface area. 
The islands can behave as an equilibrium morphological 
phasefl^j with narrow size distributions around specific 
sizes S that often do not show appreciable coarsening. 
Beyond being a recently appreciated crystal growth phe- 
nomenon, interest in the topic stems from the technologi- 
cal importance of growing strained semiconductor layers. 
Some combinations, such as InAs on GaAs, exhibit ap- 
preciable quantum confinement of electrons and holes, so 
the work has been further motivated by the study and 
possible application of island optical and electrical prop- 
erties. 

The relative contributions of growth kinetics and equi- 
librium thermodynamics to island formation have not 
been settled. Some hypotheses for the growth kinetics 
selecting the island size include that a strain-related dif- 
fusion barrier prevents atoms from moving onto larger 
islands |||7]] or that a nucleation barrier for additional 
lattice-planes on island facets limits growth J8|. An 
equilibrium-based idea for explaining the island size and 
size distribution is that nearly uniform one-monolayer- 
high precursor islands form at equilibrium and then col- 
lapse and roll up to give the observed islands 0. In 
another approach to explain the elongated islands found 
in some instances, for example Ge islands on Si p], one 
model fixes the island height and minimizes the sum of 
strain and surface energy as a function of size and as- 
pect ratio to give symmetric islands at an equilibrium 



size and elongated islands when the islands grow beyond 
the preferred size M] . 

We report on GaAs islands grown on GaP, which can 
be viewed as the III-V semiconductor analog of Ge on Si, 
since the respective lattice constants are very close to one 
another. The distinguishing behavior of this new system 
is that very long islands assemble into a dense, stable, 
corrugated phase with specific orientation and specific 
lateral period. Atomic force micrographs for samples 
with different GaAs deposition thicknesses, but other- 
wise made with the same growth temperatures, pressures 
and rates, show that the corrugations have appeared at 
a thickness near three monolayers and still dominate the 
morphology at approximately seven monolayers of de- 
position. The 13 nm corrugation period is in the range 
desired for semiconductor quantum wires with strong car- 
rier confinement to one dimension. Electron diffraction 
from the corrugated phase shows GaAs crystals com- 
pressed to the GaP lattice constant in the long direction, 
having the GaAs lattice constant perpendicular to the 
wire in the substrate plane, and being extended beyond 
the GaAs lattice constant perpendicular to the substrate. 
The present wire phase is not found to be predicted by 
an elasticity model for strained corrugated surfaces [[13| , 
and we argue that the out-of-equilibrium requirement for 
the long islands of Ref. jio) are not met under the present 
conditions. After considering the several reported orien- 
tations of other elongated islands, we conclude that the 
intrinsic surface stress anisotropy does not select island 
orientation but that inequivalcnt facet energies do. We 
modify the model of Ref. to account for the wire 
phase orientation by including inequivalent island facet 
surface energies. 

The experiments consisted of seven samples with GaAs 
depositions of 0, 0.84, 1.7, 3.3, 5.0, 6.7, and 10 mono- 
layers made in a chemical beam epitaxy machine of our 
own design. In this growth technique, molecular beams 
of chemical precursors impinge on a heated substrate 
for reaction and epitaxial growth. For both GaAs and 
GaP, the growth rate is determined by the arrival of 
Ga atoms while a slight overpressure of the respective 
As or P molecular beam is maintained. The precur- 
sors were triethylgallium (TEGa) and thermally cracked 
feriiary-butylarsine (TBAs) and tertmr?/-butylphosphinc 
(TBP). For all samples a 180 nm GaP buffer layer was 
first grown on flat (±0.2°) S-doped (001) GaP substrates 
at 640° C. The GaAs growth conditions were chosen to 
emphasize equilibrium behavior; the temperature was 
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580°C and the rate was 0.056 monolayers/s. Reflection 
high-energy electron diffraction (RHEED) intensity pat- 
terns were recorded during each growth |Q. Atomic 
force microscopy (AFM) measurements in tapping mode 
were done outside the vacuum system. 

FIG. 1. AFM micrographs for different GaAs deposition 
amounts, a. Epitaxial GaP surface (0 monolayers of GaAs). 
b. 0.84 GaAs monolayers showing small elongated islands, c. 
Corrugated surface after 5.0 GaAs monolayers, d. Irregularly 
shaped islands after 10 GaAs monolayers. 

The AFM images shown in Fig. 1 follow the sur- 
face morphology with increasing deposition, displaying 
data from the 0, 0.84, 5.0, and 10 monolayer samples. 
Monolayer-high steps are apparent on the GaP surface 
of Fig. 1(a), which was not exposed to the As molecular 
beam. At 0.84 monolayers of GaAs, Fig. 1(b), islands 
have formed with heights of up to 6 monolayers, or 2 
nm. The islands do not exhibit a sharply defined size or 
shape, though there is a clear tendency for elongation in 
the [110] direction. The widths in the [110] direction are 
20 to 27 nm, and many islands have in-plane aspect ratios 
close to 4:1. The density of islands is lxlO 10 cm" 2 . The 
density and length of the islands increase with deposi- 
tion until condensing into the corrugated phase shown in 
Fig. 1(c) for the 5.0 monolayer deposition sample. The 
corrugation spacing varies from 10 to 15 nm. Many of 
the GaAs wire structures can be seen to extend for more 
than 0.25 /im. 




Monolayers of GaAs 

FIG. 2. The first autocorrelation peak of the surface fea- 
ture heights as a function of deposited amount of GaAs. The 
peak is a measure of the separation of the nearest neighbor 
corrugations. The inset shows the autocorrelation from the 
5.0 monolayer sample taken perpendicular to the wires. 

The measured heights of the corrugations are close 
to 2 nm. If the actual island heights are greater, the 
close lateral spacing could be preventing the AFM tip 
from reaching the substrate between them. Though not 



shown in Fig. 1, the corrugations have appeared and al- 
ready cover approximately 90 % of the surface for the 3.3 
monolayer deposition. The corrugations similarly dom- 
inate the morphology of the 6.7 monolayer sample, also 
not shown here, though irregular compact islands have 
begun to form. Larger, irregular compact islands with 
heights of up to 13 nm cover the 10 monolayer sample of 
Fig. 1(d). 

Figure 2 summarizes autocorrelations of AFM feature 
heights taken for all of the samples. Autocorrelation 
was chosen to facilitate the comparison of island sepa- 
rations in non-periodic island morphologies with those 
in the approximately periodic corrugations. Height au- 
tocorrelations were calculated for each row (column) of 
data making up the image and then summed to give ID 
average autocorrelations parallel (perpendicular) to the 
wires. As an example, the result for the corrugated sur- 
face of Fig. 1(c) perpendicular to the wires is given in 
the inset of Fig. 2. The approximate periodicity in the 
inset reflects the approximate periodicity of the corruga- 
tions. The first peak in the autocorrelation is a measure 
of the average nearest neighbor separation of the islands 
parallel and perpendicular to the wires and is plotted in 
Fig. 2 for all samples. The 40 nm and 70 nm peaks of the 
0.84 monolayer sample reflect the 10 10 cm~ 2 density of 
isolated islands, where the asymmetric values are due to 
the elongated island shapes. The [110] island-island sep- 
aration decreases with deposition until the appearance 
of the corrugations near 3 monolayers, where it reaches 
13 nm and remains constant until around 7 monolayers. 
Throughout the corrugated phase the peak in the [110] 
direction decreases as the number of compact irregulari- 
ties along the wires increases. At 10 monolayers the first 
peaks are at 60 nm for both the [110] and the [110] di- 
rections, reflecting the statistical isotropy of both shape 
and positional correlations of the large irregular islands. 
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GaAs thickness (monolayers) 
FIG. 3. Diagonal strain tensor components of the GaAs 
material measured with RHEED and plotted as a function 
of GaAs layer thickness for three different crystal directions. 
The inset depicts the strain component orientations with re- 
spect to the islands. 
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The RHEED patterns recorded during growth give a 
measure in reciprocal space of the evolving surface mor- 
phology H]. When diffraction was recorded with the 
electron beam perpendicular to the wires, the diffrac- 
tion pattern changed continuously from the reflection- 
diffraction rods of a 2-fold reconstructed surface to the 
transmission-diffraction spots of the GaAs islands. With 
the beam parallel to the wires, a similar continuous tran- 
sition was recorded from an initial 4-fold pattern, though 
chevron-shaped features at the transmission spots were 
also seen. The chevron angle of approximately 35° would 
correspond to island facets along the wires of 114 orien- 
tation. 

The position of the diffraction features depends on the 
island strain. Estimates of the diagonal strain tensor 
components of the islands are plotted versus deposition 
thickness in Fig. 3, where GaP rod separations were 
used as references. For strains parallel to the surface, 
epsilon xx perpendicular to the wires and e yv along the 
wires, the (10) to (10) rod separation was measured at 
positions corresponding to 3D diffraction. Transmission 
spot separations perpendicular to the surface gave the 
strain perpendicular to the surface, e zz . Along the wires, 
in the [110] direction, the islands stay strained to close 
to the GaP substrate throughout the recorded 10 mono- 
layers of deposition. In the [110] direction, perpendic- 
ular to the wire axis, the lattice constant begins to re- 
lax to that of GaAs after approximately one monolayer, 
once the islands have formed and are contributing a sig- 
nificant amount of 3D diffraction intensity. By 3 or 4 
monolayers, the sides of the islands have bulged out to 
close to the GaAs lattice constant. The [001] tensile is- 
land strain perpendicular to the surface in the corrugated 
phase is ~0.06. Taken together, the strain measurements 
in the three perpendicular directions indicate that the 
long islands in the corrugated phase are elastically de- 
formed, implying no crystal dislocations. The present 
estimates treat the diffraction feature positions as giving 
a measure of uniformly strained islands. A more accurate 
treatment would model the nonuniform strain in the is- 
lands, the non-uniform penetration of the electron beam 
into the islands, and dynamical scattering processes dur- 
ing diffraction. The Arrows in Fig. 3 points out oscil- 
lations as function of monolayer completion, wich have 
been reported previously for the InAs/GaAs system by 

y. More detailed RHEED analysis can be found in Ref. 

The period found for the dense GaAs wire phase re- 
ported here can be compared with a continuum elas- 
ticity model that relates lattice mismatch with surface 
corrugation period (Tq] . In that model, a local mini- 
mum was found in the Helmholtz free energy, given by 
Et = E f ace t s + E ec i g es + E elastic- Here E f ace ts is the free 
energy of the facets, Edges is an energy associated with 
the facet edges and E e iastic is an elastic strain energy due 
to edges and lattice mismatch. A local energy minimum 
corresponding to the corrugated surface was found only 
to exist for a lattice mismatch that is less than a critical 



value determined by the corrugation size and material 
parameters. The critical mismatch is 
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where t xx is the facet intrinsic surface stress perpendicu- 
lar to the edges, L is the period of the faceting, and c is a 
constant dependent on the elastic moduli and geometry. 
The present GaAs corrugated surface has Lq — 13 nm. 
Here we use the same approximate value of t xx = 100 
meV/ A as in Ref. |l3| and find a critical lattice mis- 
match of (4 2 ) c = 2.3 x 10 -3 . This upper bound is an 
order of magnitude smaller than the present GaAs/GaP 
lattice mismatch of (^f) = 3.7 x 10~ 2 . Although the 
corrugated surfaces treated by the model have obvious 
similarities to the present GaAs surface, the large mis- 
match and short period mean that this model can not be 
directly applied. 

The shape of the long narrow islands of GaAs on 
GaP can be compared with a model for long rectan- 
gular islands. The island energy was taken as Ex — 
E facets + Eeiastic 0, and changes in size and in-plane 
aspect ratio were considered. For a fixed island height, 
an equilibrium size was found, and the shape of this op- 
timal island was a symmetric square. However, if islands 
larger than the equilibrium size were considered, then 
the optimal shape was found to be asymmetric with ap- 
proximately constant width but increasing length with 
size. It was suggested that isolated islands become over 
ripe during growth when they are too far apart to ex- 
change atoms through diffusion. When the height is not 
fixed, the square island always has the lower energy [|l4| . 
Consistent with the model requirements, the height of 
the GaAs islands remains approximately constant when 
they form at low coverage and throughout the corrugated 
phase. In accord with the model prediction, the width of 
the islands remains approximately constant while their 
length increases, apparently only limited by irregulari- 
ties. Still, the 13 nm corrugations are small enough that 
islands should easily exchange atoms through diffusion. 
For comparison, appreciable diffusive mass transport oc- 
curs over distances greater than 100 nm during GaAs 
growth on lithographically patterned substrates. At our 
slow deposition rate and warm temperature, the islands 
should easily exchange atoms and the model would pre- 
dict symmetric islands. Thus, the corrugated surface rep- 
resents an optimal equilibrium shape in itself without any 
over ripening. The result of Ref. |n]] does not account 
for the shape of these islands because the predicted shape 
would be symmetric for equilibrium conditions. 

The orientation of the wires might be determined by 
several contributions. One possibility is the strongly 
anisotropic surface stress which is aligned with the dimers 
of the reconstructed (001) surface. On Si, where the 
dimer orientation rotates by 90° with each atomic sur- 
face step, the elastic interaction is sufficient to sponta- 
neously rearrange a flat surface into alternating domains 
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|L7|. Here we make orientation comparisons with other 
elongated island systems. The analogous Ge islands on 
Si are 2 to 4 nm high, 25 nm wide, and can be over 100 
nm long ||. Despite the similar dimensions, the Ge is- 
lands grow with two orientations, along both [010] and 
[100], which are 45° to the present [110] oriented GaAs 
islands. In another III-V system, InSb islands on InP 
(10.4% mismatch) show a shape transition with increas- 
ing size to elongation along [110], where the islands are 
10 nm high, 30 nm wide, and 90 nm long pq| . As a 
closer mismatch analog (3.2%) having the same column- 
V elements, InAs islands formed on InP substrates show 
a shape transition with larger symmetric islands trans- 
forming to smaller elongated islands. The InAs islands 
on InP are 3 to 4 nm high, 30 nm wide, over 100 nm 
long, and aligned in the [110] direction For both the 
III-V and IV semiconductors, the dimers are parallel to 
(110) directions. The III-V materials typically have V- 
tcrminated surfaces with dimer axes parallel to the [110] 
direction. Because islands with many orientations with 
respect to the intrinsic surface stress field are reported, 
we conclude that the stress anisotropy does not play the 
determining role in island orientation. 

The difference in facet surface energies may be picking 
the island orientation. The Ge island facet energies is- 
lands will reflect the cubic crystal symmetry, and consis- 
tent with this two perpendicular orientations are found. 
The facet energies of the polar III-V compounds will re- 
flect the lack of inversion symmetry, with the largest dif- 
ference between III- and V-terminated {111} surfaces. 
While the model of Ref. jl0| does not account for long 
equilibrium shapes, the model can be modified to account 
for inequivalent orientations by introducing two facet sur- 
face energies, Ea and Eb- The the total energy divided 
by the island volume is then 

V sine/ \ t s ) 

- 2c/l ( S " lln i +t_lln ^)' (2) 

with E s the top facet surface energy; c a constant con- 
taining clastic moduli; s, t, and h the two island side 
lengths and the height; and <f> a constant containing the 
facet angle 8. The inequivalent facet surface energies 
that appear in the first term serve to provide the ori- 
entational bias. As expected, the longer island side will 
have the lower surface energy. 

To conclude, a dense GaAs wire phase was found 
to form on GaP. Diffraction patterns show elastically 
strained islands with diffraction features consistent with 
reflection from island facets. The long GaAs islands are 
similar to Ge islands found on Si. Though a crucial dif- 
ference between the two systems may be that the GaAs 
wires grow only in one orientation, allowing the conden- 
sation into a dense ordered phase. The orientation of 
elongated islands was deduced to be due to differences in 
facet energies rather than to surface stress anisotropies. 



After comparison with available theories, open questions 
include: What picks the size and what gives the long 
shapes? 
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